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Abstract 



We consider the hadronic radiation patterns for the generic process of bb + 
2 forward jet production at the LHC, where the (centrally produced) bb originate 
' either from a Higgs, a Z or from standard QCD production processes. A numer- 

ical technique for evaluating the radiation patterns for non-trivial final states is 
^ ■ introduced and shown to agree with the standard analytic results for more simple 

' processes. Significant differences between the radiation patterns for the Higgs sig- 

nal and the background processes are observed and quantified. This suggests that 
hadronic radiation patterns could be used as an additional diagnostic tool in Higgs 
searches in this channel at the LHC. 



1 Introduction 



The distribution of soft hadrons or jets accompanying energetic final-state particles in 
hard scattering processes is governed by the underlying colour dynamics at short distances 
[H 121 01 HI • The soft hadrons paint the colour portrait of the parton hard scattering, and 
can therefore act as a 'partonometer' [Tll21EllllElinilZllSliniliniinilI21. Since signal 
and background processes at hadron colliders can have very different colour structures 
(compare for example the s-channel colour singlet process gg ^ H ^ bb with the colour 
octet process qq ^ g* bb), the distribution of accompanying soft hadronic radiation 
in the events can provide a useful additional diagnostic tool for identifying new physics 
processes. 

Examples that have been studied in the literature in this way include Higgs (THj, Z' 
jO] and leptoquark [T3j production. In each case the new particle production process 
was shown to have its own particular colour footprint, distinctively different from the 
corresponding background process. 

Quite remarkably, because of the property of Local Parton Hadron Duality (see for 
example Refs. [21 El CSj) the distribution of soft hadrons can be well described by the 
amplitudes for producing a single additional soft gluon. The angular distribution of soft 
particles typically takes the form of an 'antenna pattern' multiplying the leading-order 
hard scattering matrix element squared. Confirmation of the validity of this approach 
comes from studies of the production of soft hadrons and jets accompanying large Et 
jet and ly+jet production by the CDF [J^ and DO collaborations ^Zj at the Fermilab 
Tevatron. 

One of the most important physics goals of the CERN LHC pp collider is the discovery 
of the Higgs boson. Many scenarios, corresponding to different production and decay 
channels, have been investigated, see for example the studies reported in Refs. [IHl dl 1201 • 
In a recent paper we have studied Higgs production via vector boson fusion at the 
LHC, qq — > qqH, where the colour-singlet nature of the V*V* — > H production process 
naturally gives rise to rapidity gaps between the centrally produced Higgs and the forward 
jets^. The most delicate issue in calculating the cross section for processes with these 
rapidity gaps concerns the soft survival factor 5*^. This non-universal factor has been 
calculated in a number of models for various rapidity gap processes, see for example [221 
and references therein. Although there is reasonable agreement between these model 
^Another topical example concerns central production of new heavy objects (Higgs, SUSY particles 

etc.) at hadron colliders in events with double rapidity gaps, p'p' X + M + Y (where + indicates a 
rapidity gap), which are caused by the pomeron exchanges in the i-channel. For a recent discussion and 
a list of references see |23| . 
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expectations, it is always difficult to guarantee the precision of predictions which rely 
on soft physics. However, in ^] we argued that the calculations of 5*^ can be checked 
experimentally by computing and measuring the event rate for a suitable calibrating 
process, for example the production of a Z boson with the same rapidity gap and jet 
configuration as for the (comparatively light) Higgs, see also [2H]- 

In this paper we adopt a different approach to the same problem. Rather than con- 
sidering the case where the emission of soft hadrons between the jets and central Higgs is 
suppressed (rapidity gaps) , we instead discuss the inclusive distribution using the antenna 
pattern approach to quantify the relative amounts of soft hadron emission in the signal 
and background events. In other words, we quantify how 'quiet' the signal events are 
compared to the otherwise irreducible background events. 

Thus we have in mind the following type of scenario. Suppose an invariant mass peak 
is observed in a sample of (tagged) bb events in which there are energetic forward jets, 
typical of the vector-boson production process. If such events do indeed correspond to 
Higgs production, then the distribution of accompanying soft radiation in the event — 
which we take to mean the angular distribution of hadrons or 'minijets' with energies of 
at most a few GeV, well separated from the beam and final-state energetic jet directions 
— will look very different from that expected in background QCD production of 66 + 2 jet 
events. Again, the analogous process of bb) + 2 jet production can be used to 

calibrate the analysis, since these events are, as we shall see, also generally quieter than 
the QCD backgrounds. 

Thus in this study we will consider the hadronic radiation patterns for the generic 
process of 66 + 2 forward jet production, where the (central) 66 originate either from a 
Higgs, a Z or from standard QCD production processes. We will chose configurations (i.e. 
cuts on the rapidities and tranvserse momenta of the final-state particles) that maximise 
the Higgs signal to background ratio, see j2I]. 

The paper is organised as follows. In the following section we consider the antenna 
patterns for Higgs and Z production accompanied by two forward jets. We show that 
for these colour — singlet production processes, fairly simple analytic expressions can be 
derived. However this is not the case for the more complicated QCD background processes. 
In Section 3 we show that the radiation patterns for these can be calculated using a 
more general numerical technique, which indeed can be applied to arbitrarily complicated 
processes. Section 4 summarises our results and presents our conclusions. 
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2 Hadronic Antenna Patterns for Higgs and Z + 2 Jet 
Production 

2.1 Higgs and Electroweak Z Production 

The signal process we are interested in is Higgs production via vector boson fusion, shown 
in Fig. with subsequent decay of the Higgs to 66. Furthermore, we restrict our consid- 




Figure 1: Higgs production via electroweak vector boson fusion. 



erations to the case where the outgoing quark jets are forward in rapidity and the Higgs 
decay products are central in the detector. Throughout this paper, we work in the zero 
width approximation for the Higgs and Z. As vector boson fusion involves no colour flow 
in the t-channel, the radiation pattern is simply that of the 2 — 2 process gg' qq\ with 
an additional colour disconnected hh. These were calculated in jHj. Note also that we work 
with massless quarks. The radiation pattern is defined as the ratio of the 2 — n + 1 and 
2 — i> n matrix elements using the soft gluon approximation in the former. The dependence 
on the soft gluon momentum k then enters via the eikonal factors ('antennae') 510] 

r--i_ Pi-Pj 



{Pi ■ k)iPj ■ k)' 

For the signal q{pi) q'{p2) — ^ liPs) <i'ijPi) + gik) we have 



:2.i) 



\M? = '^{ '^^] 2C,m + [24]) (2.2) 



with S = (pi + P2)^ t = {Pl- P3)^ U={pi- P4)2. 

We then normalise this by the matrix element for the leading order process q{pi) q'{p2) - 

(i{p3)q'{p4)- 

ElA^I^ - (f!±^) . ,2.3, 
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Note that in this particular case, the 2^3 matrix element in the soft gluon limit factorises 
into the form (2 — > n matrix element) x (antenna factor). This feature is not universal, 
being restricted to only very simple cases such as this. The antenna pattern is then 

7J(,,' ^ = _ ^C'HllS] + [24]). (2.4) 

As we are working in the zero width approximation^ we can include the decay of the Higgs 
into (massless) 66 by simply adding the antenna for this colour disconnected part. The 
hadronic radiation pattern for q{pi) q'{p2) — q{p^) q'{pi) H ; H ^ Kp^) Kps) is then 

n{H) = 2Cf ([13] + [24] + [56]). (2.5) 

In order to visualise the pattern we must specify the kinematics and pick a relevant 
configuration for the incoming and outgoing particles. We label the four-momenta by 

a{pi) + h{p2) c{p3) + d{p4) + ■■■ + g{k), (2.6) 

where the gluon is soft relative to the other large-i^T final state partons, i.e. k ^ Et- 
We ignore the gluon momentum in the energy-momentum constraints, work in the overall 
parton centre of momentum frame, fix the Higgs to be at rest in that frame and its decay 
products at (?7b, 0b) = (0, 7r/2) and (0,37r/2). With the notation p'^ = {E,px,Py,Pz), the 
momenta are then 

Pi = {■mH/'2 + ETCOshrijet,0,0,mH/'2 + ETCOshrjjet), 

P2 = (m/^/2 -F£;t cosh ryjet,0, 0,-771^/2 --Et cosh r^jet), 

p^ = (Et cosh rijet,0,ET,ET sinh r^jet ) , 

p1 = {Et cosh?7jct, 0, -Et, -Et sinh rj^et), 

p>i, = (m^,, 0,0,0), 

p^ = (m^^/2,mH/2,0,0), 

pI = (m^^/2,-m^/2,0,0), 

— {kT cosh r)g, kT sin 0g, /cy cos (pg, kT sinh rjg). (2.7) 

This is the appropriate form for studying the angular distribution of the soft gluon, 



^Actually our analysis is formally correct provided that Th ^ Eg where Eg is the typical soft 
gluon/hadron energy, i.e. the Higgs lives long enough to prevent any interference between gluon emission 
before and after the Higgs decays. In any case, such interference would occur only in higher orders in Ug 
and is colour suppressed. 
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parametrised by rjg and (pg. Using the kinematics of Eq. ()2.7|1 with Eq. ()2.5|1 gives 
2Cf i coshr^jet - sinhT^jet 



+ 



+ 



k'^ \ (cosh ?7g — sinh ?7g) (cosh T^jet cosh r]g — cos 0g — sinh r^jet sinh rjg 

cosh ?7jct - sinh rj^^t 

(cosh rig + sinh r/g) (cosh r/jet cosh r/g + cos (pg + sinh r^jet sinh T^g) 

2 



(2 8) 

(coshr^g — sin(0g + 7r))(cosh?7g + sin(0g + vr)) J 

Note that the result is independent of Et and ttih and that colhnear singularities are 
situated at (%,0g) = (?7jet, tt) , (— r^jct, 0) , (0, 7r/2) and (0, 37r/2). As an illustration, Fig- 
ure |21 shows k^TZ{H) with r^jet = 3.5. One can clearly see that a colour connection exists 




Figure 2: Antenna pattern for qq' — > qq'H ; H ^ bb. Here rjjet = 3.5. 



Figure 3: The point {r]g,(pg) = (0, 7r/2) in lZ{qq' qq') as one varies r/jet- As the jets move 
apart, the antenna falls to zero. 

between the initial state parton pi and final state jet p3, similarly with p2 and J94, and 
also between the 6-quark jets. The antenna pattern is small between the jets and the 6's 
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as there is no colour connection between these — this is the 'rapidity gap' phenomenon. 
The emission of soft gluons in the rapidity gaps decreases as the gap widens. This is 
illustrated in the case without the 6-quark antenna (Fig. EJ, which shows the antenna 
pattern at (?7g, 0g) = (0, 7r/2) as a function of ?7jet- 

qA q,q z q q 




q,q q,q q,q q,q q z 

(a) (b) (c) 

Figure 4: Electroweak Z production. 

Next we consider the analogous electroweak Z production process (Fig. |3)), which can 
in principle be used to calibrate the Higgs production process. In this case the variety of 
diagrams at leading order means that there is no exact eikonal factorisation. However in 
the kinematic limit we are interested in — forward jets and central Z production — the 
dominant amplitude is again the one involving t-channel W exchange, i.e. WW Z, and 
the antenna pattern is trivially identical to that for Higgs production. We will prove this 
result, and consider its implications, when we discuss how to calculate antenna patterns 
numerically below. 

2.2 QCD Z Production 

In practice, Z + 2 jet production can also occur by O^a^aw) QCD production involving 
t-channel gluon exchange, see Fig. Because of the different colour structure of such 
diagrams we would expect a very different antenna pattern. 

Once again there is no exact factorisation of an overall soft gluon form factor and 
therefore no simple expression for the radiation pattern. However, as for electroweak Z 
production the factorisation is restored in the forward jet - central Z limit, in which 
case the antenna pattern is identical to that for the QCD 0(a|) qq' —>■ qq' production 
process |9J, i.e. 

n{QCD Z) ^ 2Cf{M + [23]) + ^[12; 34] + 2Cf[56], (2.9) 
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where 

[tj; kl] ^ 2[zj] + 2[kl] - [zk] - [zl] - [jk] - [jl]. (2.10) 

Substituting the kinematics of Eq. ()2.7|) and plotting the resulting analytic expression 
with rijet = 3.5, one obtains Figure IHl Before commenting on the differences, we note that 
both Figures El and IHl exhibit the same limiting behaviour 

lim k^n{H,QCD Z) = ACf, (2.11) 

as a consequence of both processes having initial state quarks^. They are also identical 
as one approaches the colhnear singularities corresponding to the final state 6-jets: 

hm e^n{H,QCD Z)^ACf—-^ — -— . (2.12) 

(%=9ig)^('?jct,9ijct) cosh (?7g - ?7jet) - COS^(0g - 0jet) 

The difference in the colour flow shows up in the region between the two final state forward 
quark jets, as expected. Taking the ratio of the two patterns makes this difference plain 
(Fig. [7j). The maximum difference occurs at (?7g,0g) = (—4.4,0) and (r/g, 0g) = (4.4, tt) 

■^Of course qg — > Zqg also contributes to Z + 2 jet production, and this will have a different colour 
structure from qq — > Zqq. For purposes of comparison with the Higgs case, we only consider quark 
induced production in this section. 
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when the ratio attains the value 2.3. This shows the colour connection between the initial 
state (at rjg = ±00) and the forward jets in the Higgs production case that is suppressed 
by a factor 0{-^) in the QCD Z-production case. Another interesting phase space point 
is at 0g) = (0,0), i.e. the central region transverse to the bb axis. Here the radiation 
pattern increases by a factor of three going from Higgs to QCD Z production, indicating 
the presence of an additional underlying colour connection in the latter case. 



fe|7^(Jf) 2 

kl.n{QCD,Z) 




Figure 7: Ratio of Higgs to QCD Z production antenna patterns. The ratio is unity at the 
position of the jets (indicated by arrows) and in the forward and backward limits 
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3 Numerical Hadronic Antenna Patterns 



An important (and dominant) background to the processes considered in the previous 
section comes from QCD 0{a^) 66 + 2 jet production when M^^ ~ Mh,z'^- Some sam- 
ple diagrams are shown in Fig. |H| There is clearly no unique and simple colour flow 
associated with these diagrams, and hence no compact analytic antenna pattern can be 
derived. This is an example of a situation where there is no factorisation of the form 
(2 — > 4 matrix element) x (antenna factor). However we can instead use a purely numeri- 

Q,Q Q,q 9 9 




q,q q,q 9 9 

Figure 8: Examples of QCD dijet plus bb production diagrams. 

cal method in which we compare the values of the 2 ^ n and 2 — > n + 1 matrix elements at 
each point in phase space, their ratio in the soft gluon limit deflning the antenna pattern. 
In order to verify that this methodology works, and in particular to establish how soft the 
gluon has to be before the limiting pattern is reached to some level of precision, we flrst 
make a numerical evaluation of the analytic radiation patterns discussed in the previous 
section. 

3.1 Comparison of Numerical and Analytic Antenna Patterns 
for Signal Processes 

Unlike the analytic case, where we can simply ignore the momentum of the soft gluon 
in assigning a kinematic conflguration that respects momentum conservation, we must 
account for the numerically flnite gluon momentum in evaluating the matrix elements. 
Thus there is a degree of arbitrariness introduced. We choose to assign the momenta such 
that the central boson or bb system cancels the 3-momentum of the soft gluon. In other 
words 

= y^lH + k\ -k) . (3.1) 

*We are not discussing here the background caused by a possible misidentification of the gluons as b 
jets. For a recent treatment of this see |27| . 
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Therefore the value of the antenna pattern depends on the specific kx that we choose for 
the soft gluon, but in such a way that k^TZ tends to a finite hmit as fc^ —>■ 0. Figure IHl 
illustrates this by taking the ratio of the numerical qq' —>■ qq'H antenna pattern with the 
analytic qq' qq' antenna pattern for kxg = I GeV. The ratio is close to unity, except 
when the gluon rapidity is very large. In this region the 'soft' gluon carries a significant 
amount of energy and begins to distort the overall kinematics. For numerical purposes 
only, as a formal check that this effect is under control, we can set kxg to be sufficiently 
(and artificially) small to make sure the analytic result is recovered everywhere. Thus 
Fig. El shows the same ratio for kxg = 10~^ GeV - no deviation from unity is now 
discernible. Note that we will always use kxg = 1 GeV in making predictions for the 
antenna patterns using the numerical treatment. Since our ultimate aim is to compare two 
numerically generated antenna patterns in signal to background studies, the discrepancies 
at large gluon rapidity visible in Fig. IHl will exactly cancel in the comparison. 



njqq'^qq'H) 
niqq'^qq') 




Figure 9: Ratio of numerical qq' — > qq'H to analytic qq' qq' antenna patterns with |r;jet| = 2 
and kxg = 1 GeV. 

As already pointed out, the antenna pattern for the full electroweak qq' —>■ qq' Z 
process is not given by the simple analytic approximation, except when the jets are far 
forward. We can now illustrate this using the numerical method. Thus Figs. ^2 and IT^ 
show the ratio of the numerical electroweak qq' qq' Z antenna pattern with the analytic 
electroweak qq' qq' antenna pattern for the choice of l^^jetl = 2 and |?7jet| = 4 with 
kx = 10"^ GeV. In the former case, the agreement with the analytic antenna pattern is 
only at the 10% level. The discrepancy is due to the contribution of the Z-sstrahlung 
diagrams (Fig. |31d) in the numerical case. However, as one forces the quark jets to be 
more forward the discrepancy decreases. Therefore, as long as we require the jets to be 
forward (i.e. \t\ <^ y/s), the analytic approximation is valid. 
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I.ln 

1.0S- 

n(qq'^qq'ZQCD) ^,m- 



H{qq'^qq') 




Figure 13: Ratio of numerical QCD qq' Qq'Z to analytic qq' — > qq' antenna patterns with 
\r]jet\ = 2 and kxg = 10~^ GeV. 
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Figures ITT?1 and IT^ show the same quahtative effect in the QCD mediated Z production 
case. The deviation from our approximation that \t\ ^ \fs is noticeably less than in the 
electroweak case. The reason for this is that in the electroweak case we are kinematically 
disturbing a delicate interplay between the numerator and the denominator in the term 
describing the colour connection between pi and 



[13] 



(3.2) 



(pi ■ k) {p3 ■ k) 

In particular, due to the smallness of the numerator, this contribution is strongly sup- 
pressed for the radiation outside the narrow cones around the directions of the incoming 
and outgoing partons. Contrast this with the QCD Z production case where the dominant 
colour connection is between pi and p4 

Pi - Pi 



[14] 



(3.3) 



{Pi ■ k) {p4 ■ k) 

Here the numerator is not small. This cancellation is therefore more stable and our 
kinematic disturbance has less effect. 



1.1 



njqq'^qq'ZQCD) i,oB 
■Riqq'^qq') 



1,M 



1,02 




Figure 14: Ratio of numerical QCD qq' qq'Z to analytic qq' qq' antenna patterns with 
|?7jet| = 4 and kxg = 10"^ GeV. 



3.2 Numerical Antenna Patterns for Background Processes 

Fi gure ITHl shows the numerical antenna pattern for the QCD mediated process qq' —>■ qq'bb. 
We will again focus mainly on the background process with initial state quarks, to allow 
comparison with the signal processes. In any case, the typical of the parton-level 
process is typically several TeV at the LHC^, so we are working at high x and quark 

^For example, from Eq. \/i ~ + 2£;t cosh r/jot =i 2.8 TeV for Et ^ 50 GeV tuh = 120 GeV 

and ryjct = 4. 
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initiated processes will dominate. Therefore the antenna patterns for the signal and 
background processes become identical near the beam and final state 6-quark directions, 
being dominated by the (universal) collinear singularity for emission off quark lines. 

Figure El shows the radiation pattern for the background QCD process qg —>■ qgbb 
with |77jet| = 4 and kxg = 1 GeV. As expected, the pattern is much more complicated than 
that for the signal H or Z production processes. Colour strings can now connect many 
more pairs of initial and final state particles, and the overall level of radiation is higher as 
a result. However in the directions of the incoming and outgoing partons, the distribution 
of soft radiation is the same as that for the signal processes. Thus, in particular, the 
distribution approaches iCp for large positive rjg, cf. Eq. 1)2.111) . 

For completeness, we show in Figs.^l-Elthe corresponding antenna patterns for the 
other QCD 2^2 + (bb) processes. The most obvious differences are in the size of the 
distributions near the incoming and outgoing partons, where the limiting ACp behaviour 
for emission off quarks is replaced by 40^ for emission off gluons. 

The interesting quantities are of course the differences between the signals and back- 
grounds. Figure [201 shows the ratio of numerical qq' —>■ qq'H; H —>■ bb to numerical 
qq' — >■ qq'bb antenna patterns, for the same typical kinematic configuration as before, i.e. 
I^jetl = 4 and kjng = 1 GeV. We see that the ratio (i) falls to near zero between the central 
and forward particles (rapidity gap effect), (ii) is larger than one between the final-state 
bb pair, (iii) is larger than one between the forward jets and the beam (the [13] and [24] 
connection in the signal), and (iv) approaches unity in the forward/backward directions 
and at the locations (marked as arrows) of the incoming and outgoing particles. Over the 
whole [f], 0) plot, the ratio varies in size from a minimum of 0.03 to a maximum of 2.3, 
i.e. a factor of 70. The corresponding ratio of the antenna patterns for the electroweak 
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Figure 16: Numerical antenna pattern for qg qgbb with jryjetl = 4 and kxg = 1 GeV. 




Figure 17: Numerical antenna pattern for gg ggbb with \r]jet\ = 4 and k^g = 1 GeV. 




Figure 18: Numerical antenna pattern for qq ggbb with jryjetl = 4 and kxg = 1 GeV. 
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Figure 19: Numerical antenna pattern for gg qqbb with |7/jet| = 4 and kxg = 1 GeV. 

Z production and QCD background is of course very similar. 

We next consider (Fig. I?H) the ratio of the QCD Z-production and background 
qq' — s> qq'bb antenna patterns. There is much less structure here than there was in the 
corresponding Higgs case - note in particular that the rapidity gap dip between the for- 
ward and central particles is absent^, and indeed that the ratio is close to one everywhere 
except near the central b jets. In the Z production case, there is always a colour string 
connecting the b and the b, and this results in the ratio increasing to a maximum of about 
1.5 between these two particles. This value has a weak dependence on the rapidities of 
the forward jets. Figure shows the slice through Fig. |^ at ?7g = as l^^jetl is varied from 
1 to 8. The ratio is always one at 0jet = tt and 37r/2, the location of the b and b. 

4 Conclusions 

Hadronic radiation patterns can provide a useful additional tool enabling us to improve the 

separation of Higgs production from the conventional QCD-induced backgrounds. In this 

paper we have focused on the vector boson fusion mechanism of Higgs production in the 

events with two forward tagging jets. We find that the fairly simple analytical expressions 

reflecting the coherent structure of QCD radiation off the multi-parton system (antenna 

pattern) can serve quite successfully as a qualitative guide for the more general numerical 

^Note that by imposing tlie rapidity gap requirement to isolate the centrally produced system from 
the proton remnants, we would automatically cut off the colour connection between this system and 
the forward going partons. As shown in this allows us to substantially reduce the background 
contributions, though at the price of a reduction in the overall event rate (due to the notorious survival 
factors) . 
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Figure 20: Ratio of numerical qq' — > qq'H; if —> 56 to numerical qq' — > qq'bb with \r]jet\ = 4 
and kxg = 1 GeV. 
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Figure 22: Slice in r/g = of Fig. as the rapidity of the forward jets is varied. 
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calculational technique, which in turn can be apphed to a large variety of comphcated 
processes. 

The analysis presented here should be regarded as a 'first look' at the possibilities 
offered by hadronic flow patterns in searching for the Higgs in vector boson fusion. Of 
course, ultimately there is no substitute for a detailed Monte Carlo study including detec- 
tor effects. However the results presented here indicate that the effects can be potentially 
large, and therefore that more detailed studies are definitely worthwhile. 
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